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A NUMERICAL METHOD FOR CONSTRUCTING THE 
\/V HYPERBOLIC STRUCTURE OF COMPLEX HENON MAPPINGS 

o: 
o- 

{SJ ■ SUZANNE LYNCH HRUSKA 

o : 

Q Abstract. For complex parameters a, c, we consider the Henon mapping 

Ha,c : C 2 — > C 2 , given by (x,y) h-> (x 2 + c — ay,x), and its Julia set, J. In this 
paper, we describe a rigorous computer program for attempting to construct 
a cone field in the tangent bundle over J, which is preserved by DH, and a 
continuous norm in which DH (and DH~ 1 ) uniformly expands the cones (and 
their complements). We show a consequence of a successful construction is a 
rf\ • proof that H is hyperbolic on J. We give several new examples of hyperbolic 

Qmaps, produced with our computer program, Hypatia, which implements our 
methods. 

•S: 

i— — i. 1. Introduction 

£NJ ■ The Henon family, Ha, c {x,y) = (x 2 + c — ay,x), has been extensively studied 

^ j as a diffcomorphism of R 2 , with a, c real parameters. For example, Benedicks and 

Carlcson show the existence of chaotic behavior in the form of a strange attractor for 
^^ | some real Henon maps in jHj. Here we consider H ac as a diffeomorphism of C 2 , and 

iq ■ allow a, c to be complex. Foundational work on the dynamics of the complex Henon 

O ! family has been done by Bedford and Smillic (0EH3IH), Hubbard ([231211122), 

and Fornaess and Sibony f[14|). However, basic questions remain unanswered. 

A natural class of maps to study are the hyperbolic maps, since hyperbolic 

r^ ■ maps generally have nontrivial (chaotic) dynamics, but are amenable to analysis. 

"ti | A Henon mapping H a . c is hyperbolic if its Julia set, J a ,c, is a hyperbolic set for 

r^ ■ H a .c- (Hyperbolicity and Julia set are defined in Section 2.) For Henon mappings, 

hypcrbolicity implies Axiom A, which implies shadowing on J, i.e., e-pscudo orbits 

are (5-close to true orbits, and structural stablity on J, i.e., in a neighborhood in 

parameter space the dynamical behavior is of constant topological conjugacy type. 

«2j ■ Thus for a hyperbolic mapping, the dynamics on J should be able to be understood 

using combinatorial models. These properties make hyperbolic diffcomorphisms 

amenable to exploration via computers. 

Motivated by careful computer investigations, Oliva ([221) provides a combina- 
torial model of the dynamics of some Henon mappings, including for example, the 
mapping of Figure Q The proposed model presupposes that the mapping is hyper- 
bolic. Hubbard and Papadantonakis (^GO!) have more recently generated pictures 
of slices of the Henon parameter space, which attempt to sketch either the locus 
of maps with J connected, or the locus of maps with J having no interior (see 
Figure [2J|. These and other computer investigations suggest that the dynamical 
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Figure 1. The filled Julia set for a Henon mapping iJ Q;C (x, y) = 
(x 2 +c — ay, x), restricted to the unstable manifold of a saddle fixed 
point, with its natural parameterization. Here a = .3, c = —1.17, 
and the map has attracting cycles of periods one and three, which 
is impossible for quadratic polynomial maps of C 




Figure 2. A slice of Henon parameter space, the c-plane, with 
a = .3. The innermost dark set is an attempt to sketch the set of 
maps with connected Julia sets, while the exterior is an attempt 
to sketch the maps with filled Julia set having empty interior (in- 
cluding complex horseshoes). These regions have features reminis- 
cent of the parameter regions of one-dimensional polynomial maps. 
This is an intriguing parallel, suggesting that the Henon parameter 
space may be an equally rich arena of study. 
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behavior of the complex Henon family is rich and subtle. If certain of these map- 
pings could be shown to be hyperbolic, then this would serve as a first step toward 
showing mathematically that these apparent phenomena actually occur. 

However, there are very few complex Henon mappings known to be hyperbolic. 
Let us summarize what is known. First, if 

(1) |c|>2(l + |a|) 2 , 

then -ff a ,c|j is conjugate to the full 2-shift (so J a . c is a Cantor set), and the map is 
hyperbolic (compare Devaney and Nitecki ^21> Oberste-Vorth [31], Morosawa, et 
al. |28j). In this case the mapping is called a (complex) horseshoe. The exterior 
in Figure contains the set of horseshoes (among other types of maps). Second, 
Hubbard and Oberste-Vorth ((2Hj) show that if P c (x) = x 2 + c is a hyperbolic 
polynomial, then there exists an A(c) such that if 

(2) 0< \a\ <A(c), 

then H ac \j is topologically conjugate to the function induced by P c on the inverse 
limit lim^(J, P c ), hence H ac is hyperbolic. Ishii and Smillie ([21]) have worked 
to obtain explicit estimates for the constant A(c) in J2]), but these estimates are 
relatively small. 

Our broad goal is to develop computer algorithms with which we can rigorously 
describe the dynamics of any hyperbolic complex Henon mapping. In this paper 
we make a key step in that process, by developing a computer program which 
can establish if a complex Henon mapping is hyperbolic, and if so, the program 
produces explicit information about how the map is hyperbolic; in particular, it 
builds two complementary cone fields in the tangent bundle over J (the unstable 
and stable cones), and constructs a norm in which DH (DH~ l ) preserves and 
unifomly expands the unstable (stable) cones. 

Since hyperbolicity is structurally stable, a computer program with infinite re- 
sources should be able to prove hyperbolicity for any hyperbolic map. However, 
non-hyperbolicity is an unstable condition, thus a computer program cannot be 
expected to recognize when a map is not hyperbolic. (For example, P(z) = z 2 — 1.5 
is presumably, but not provably, non- hyperbolic.) 

This paper builds on the results of [TH] (cf [23CHE3E1E])- There we describe 
an algorithm, called the box chain construction, which given e > finds a compact 
neighborhood, 23, containing a 6 = <5(e)-neighborhood of J, and creates a finite 
graph, r, which models the e-dynamics of H on 23. In this paper, we need only 
know the following about T. 

Definition 1.1. Let T = T(V, £) be a directed graph, with vertex set V = {Bk}^ =1 
a finite collection of closed boxes in C 2 , having disjoint interiors, and such that the 
union of the boxes 23 = U^ =1 Bk contains J. Suppose there is a S > such that T 
contains an edge from Bk to Bj if the image H(Bk) intersects a 5- neighborhood of 
Bj, i.e., 

£D{(k,j):H(B k )nX(B 3 ,5)^(t)}. 

Further, assume T is strongly connected, i.e., for each pair of vertices B k , Bj, there 
is a path in T from B^ to Bj , and vice- versa. Then we call T a box chain model of 
H on J. 
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In |18| . we describe how the box chain construction builds strongly connected 
graphs modeling every basic set of H, for example, J and any attracting periodic or- 
bits. In this paper, we let T be the strongly connected graph component containing 
J , and we mostly ignore the others. 

Our first task in this paper is to develop a discrete condition on T we call 60a; 
hyperbolicity, and show this condition implies hypcrbolicity of H: 

Theorem 1.2. Let T be a box chain model of H on J. IfT is box hyperbolic, then 
H is hyperbolic on J . 

Our definition of box hyperbolicity is inspired by our work in numerically estab- 
lishing hypcrbolicity in one complex dimension (see JHj)- The difference is that 
in one dimension, hypcrbolicity simply means expansion on J, whereas for Hcnon 
mappings, hypcrbolicity means a saddle property, i.e., expansion, contraction, and 
transversality. The notion of box hyperbolicity is made precise in Section 0J Let 
us briefly describe this property here. We begin with the cone field criterion for 
hyperbolicity. In a more general setting, Newhouse and Palis ([2111231) show that 
an /-invariant set A is hyperbolic for / iff there is a field of cones in the tangent 
bundle over A such that Df maps the cone field inside itself, and such that in some 
norm, Df uniformly expands the cones, and Df -1 uniformly expands the comple- 
ments of the cones. Moreover, the field of cones need not be continuous in x E A; 
hence the cone field criterion for hyperbolicity yields a natural way to study the 
hyperbolic structure of a diffcomorphism using a computer (i.e., discretely). Here, 
we build cones which are constant on each box vertex of I\ 

To use the cone field criterion, we must find both a field of cones preserved by 
DH and a norm in which DH expands the cones. We cannot expect that DH 
expands vectors in each cone with respect to the euclidean norm. For example, 
in the euclidean norm, DH may be small for some cones over a pseudo-cycle, but 
larger in others, so that only the total cycle multiplier is larger than one. Thus 
given a model T of a map H, we attempt to build a discretized norm on the tangent 
bundle over T, which is designed to factor out the differences in DH along cycles, 
so that in this new norm, DH is expanding on every cone. Then we show that 
hyperbolicity in the discrete norm implies hyperbolicity for some continuous norm. 

In order to test for box hyperbolicity on a given T, we have developed a com- 
puter algorithm we call the Axis Metric Algorithm, designed to either prove that a 
given r is box hyperbolic, or describe which parts of T arc obstructions to proving 
box hyperbolicity. If a T fails to be box hyperbolic, then either the map is not 
hyperbolic, or the boxes of T are too large. Thus our approach is to attempt to 
prove box hypcrbolicity on a sequence of models T(n) with decreasing box size. 
If T(n — 1) fails to be box hyperbolic, then to create T(n) we could decrease the 
size of all the boxes, or use output of the Axis Metric Algorithm on T(n — I) in 
choosing which boxes need to be decreased in order to create a T(n) more likely to 
be box hyperbolic. If some T(n) is found to be box hyperbolic, then / is hyperbolic, 
and the program terminates. Thus a "successful" run of an implementation of our 
procedure gives a mathematical proof of hypcrbolicity: if we can construct a T 
which the Axis Metric Algorithm shows to be box hyperbolic, then that mapping 
is mathematically proven to be hyperbolic, by Theorem 1 1.21 

In designing the Axis Metric Algorithm for verifying box hypcrbolicity of some 
r, we build on the one-dimensional procedure described in |19j . There we prove 
hyperbolicity of polynomial maps of C by creating a piecewise continuous (box 
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constant) metric, under which the map is expanding on a neighborhood of J. To 
move up to two dimensions, and saddle-type hyperbolicity, we break down the 
problem into one dimensional pieces, then reassemble. In particular, we build 
approximately invariant, box constant unstable and stable line fields, which will 
serve as axes for our cones. Then we use the one dimensional algorithm twice, to 
attempt to build a metric which is contracted on the stable directions, and another 
metric which is expanded on the unstable directions. These metrics and line fields 
then determine a cone field with an induced metric, which we use to test for box 
hyperbolicity. The Axis Metric Algorithm is described in detail in Section [SJ 

Finally, we have implemented our methods into a computer program called Hypa- 
tia, and used the program to prove hyperbolicity of several Henon mappings which 
were not previously known to be hyperbolic: 

Theorem 1.3. The complex Henon mappings, H a ^ c (x,y) = (x 2 + c — ay,x), with: 

(c, a) = (-0.3, 0.1), (0, -0.22), (-3, -0.25), and (1.5, 0.5), 

are hyperbolic. 

Computer pictures suggest that the first two mappings of Theorem 11.31 H ac 
with (c,a) = (— .3, .1) and (c,a) = (0, — .22), are in the main cardioid, with the 
only recurrent dynamics consist of a connected J and one attracting fixed point, 
and that the latter two mappings are horseshoes, with J ac with c = —3, a = —.25 
appearing to lie in R 2 , and J a , c , with c = 1.5, a = .5 appearing not to be contained 
in R 2 . Whether or not that is the case, each of the maps of Thcorcm ll.3l lics outside 
of the known regions defined by JQ) and J2J with the Ishii-Smillie estimates. 

All of the computations involved in proving Theorem 11.31 were run on a Sun 
Enterprise E3500 server with 4 processors, each 400MHz UltraSPARC (though the 
multiprocessor was not used) and 4 GB of RAM. 1 When computations became 
overwhelming, memory usage was the limiting factor. The C-\ — h, unix program, 
Hypatia, may be obtained from the author. 

To conclude the introduction, we sketch the organization of the paper. We 
give background on the dynamics of the Henon family in Section In Section 
we briefly discuss interval arithmetic with directed rounding, the method used to 
maintain rigor in our computer computations. In Section 0J we define box hyper- 
bolicity, and we prove Theorem 11.21 establishing that box hyperbolicity implies 
hyperbolicity In Section [5] we describe our computer procedure for verifying box 
hyperbolicity, including the Axis Metric Algorithm. Finally, in ScctionElwe provide 
some data on how we used Hypatia to prove hyperbolicity of each of the maps of 
Theorem 1 1.31 
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2. Background 



2.1. The Henon Family. Polynomial diffcomorphisms of C 2 necessarily have poly- 
nomial inverses, thus are often called polynomial automorphisms. Fricdland and 
Milnor (|15p showed that polynomial automorphisms of C 2 break down into two 
categories. Elementary automorphisms have simple dynamics, and are polynomi- 
ally conjugate to a diffcomorphism of the form (x, y) i— > (ax+b, cy+p(x)) (p polyno- 
mial, a, c ^ 0). Nonelementary automorphisms arc all conjugate to finite composi- 
tions of generalized Henon mappings, which are of the form f(x,y) = (p(x) — ay,x), 
where p{x) is a monic polynomial of degree d > 1 and a^O. 

To clarify the situation, one can define a dynamical degree of a polynomial au- 
tomorphism of C 2 . If deg(f) is the maximum of the degrees of the coordinate 
functions, the dynamical degree is 

d = d(f)= lim {deg{r)) X/n - 

n — >oo 

This degree is a conjugacy invariant. Elementary automorphisms have dynamical 
degree d = 1. A nonelementary automorphism is conjugate to some automorphism 
whose polynomial degree is equal to its dynamical degree. Without loss of gen- 
erality, we assume such / are finite compositions of generalized Henon mappings, 
rather than merely conjugate to mappings of this form. 

Thus, the quadratic, complex Henon family H ac (x,y) = (x 2 + c — ay, a;) rep- 
resents the dynamical behavior of the simplest class of nonelementary polynomial 
automorphisms; those of dynamical degree two. In this paper, we usually use the 
letter / for a polynomial diffcomorphism of C 2 with d(f) > 1, and H for a (degree 
two) Henon mapping. We state results in Section 0] for the more general /, but in 
explaining the procedure for verifying hyperbolicity in Sectional we concentrate on 



2.2. Drawing Meaningful Pictures. For a polynomial map P of C, the filled 
Julia set, K, is the set of points whose orbits arc bounded under P; the Julia set, 
J, is the topological boundary of K. For a polynomial diffcomorphism /, like H a ^ c , 
there are corresponding Julia sets: 

• K + (K~) is the set of points whose orbits are bounded under f{f~ l ) 
and K = K + n K~~ is called the filled Julia set; 

• J ± = dK^ (topological boundary) 

and J = J + fl J~ is called the Julia set. 

Filled Julia sets are the (chaotic) invariant sets which can be easily sketched 
by computer, on any two-dimensional slice. Hubbard has suggested the following 
method for drawing a dynamically significant slice of the Julia set, by parame- 
terizing an unstable manifold. This method has been implemented by Karl Pa- 
padantonakis in FractalAsm ([201 H]). Figures d E| and were generated using 
FracalAsm. 

Let / be a diffcomorphism of C 2 . If p is a periodic point of period to, and 
the eigenvalues X,fi of D p f m satisfy |A| > 1 > |/z| (or vice-versa), then p is a 
saddle periodic point. The large (small) eigenvalue is called the unstable (stable) 
eigenvalue. If p is a saddle periodic point, then the stable manifold of p is W s (p) = 
{q: d(f n (q),f n (p)) — > as n — > oo}, and the unstable manifold of p is W u (p) = 
{q: d(f~ n (q),f~ n (p)) — > as n — > oo}. If p a saddle periodic point of /, then 
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W u (p) (W s (p)) is biholomorphically equivalent to C, and on W u (p) (W s (p)) 1 f is 
conjugate to multiplication by the unstable (stable) eigenvalue of D p f. 

When \a\ ^ 1, except on the curve of equation Ac = (1 + a) 2 , the map H ac has at 
least one saddle fixed point, p, (|20|). The unstable manifold W u (p) has a natural 
paramctrization 7 : C — ► W u (p) given by 

7 (z) = lim J m (z)= lim H m [p+^-w x 

m — >oc 771 — >oo V Ai 

where Ai is the unstable eigenvalue of D p H and vi is the associated eigenvec- 
tor. This parametrization has the property that H( / y(z)) = j(\z), and any two 
parametrizations with this property differ by scaling the argument. 

Observe that since W u (p) C A" - , to get a picture of K in W u (p) we need 
only color pixels black which are guessed to be in K + . To sketch the picture, we 
approximate 7 by say g = 7100 in a region in the plane, B = {z = x + iy: a < 
x < b,c < y < d}. Then an escape threshold is chosen, like 10, and then for each 
z € B, if \\H n (g(z))\\ < 10 for all n < 100, we say g{z) G K+ and color it black. 
Otherwise, color according to which iterate ||-ff n (<7(.z))|| first surpassed 10. 

2.3. Hyperbolicity. In this subsection, let / be a diffeomorphism of a manifold 
M, and let A be a compact, /-invariant set. First we recall the standard definition 
of hyperbolicity (see [3"fip. 

Definition 2.1. A is hyperbolic for / if at each x in A, there is a splitting of the 
tangent bundle T X M = E x © E%, which varies continuously with 16A, such that: 

(1) / preserves the splitting, i.e., D x f(E°) = Ej x , and D x f(E%) = E^ x , and 

(2) Df (Df^ 1 ) expands on E u (E s ) uniformly, i.e., there exists a constant 
A > 1 and a norm ||-|| on T^M, continuous for x 6 A, for which 

\\D x f(vs u )\\ f(x) > X ||w u || x , for w u G E%, and 

ll^r^w 8 )!!^,^ > A||w-|| B ,forw-e^. 

As noted in the Introduction, Newhouse and Palis f |29H30| ) show hyperbolicity 
can be described using a cone field. To define a cone C x for each point x in A, we 
need a splitting T X M = E\ x © E^ x , and a positive real- valued function e(x) on M. 
Then define the e(x)-sector S e ^(Ei x ,E2 X ) by 

S e{x) {E lx ,E 2x ) = {(vi,v 2 ) G E lx ®E 2x : ||v 2 || < e{x) |[vi||}. 

Then C x = S e r x y Newhouse and Palis show that A is hyperbolic for / iff there is 
a field of cones {C x C T X M : a: £ A}, a constant A > 1, and a continuous norm ||-||, 
such that Df preserves the cone field, i.e., D x f{C x ) C Cf x , and such that in this 
norm, D x f (D x f~ l ) uniformly expands vectors in C x {T X M — C x ); moreover, the 
field of cones x — ► C x need not be continuous. Our computer algorithm for verifying 
hyperbolicity actually combines these two notions, as we will see in Section 0J 

Bedford and Smillie (|H]) have shown that for / a polynomial diffeomorphism 
of C 2 , with d(f) > 1, / is hyperbolic on its Julia set, J, iff / is hyperbolic on its 
chain recurrent set, 3^, iff / is hyperbolic on its nonwandcring set, f2. Thus we say 
/ is hyperbolic if any of these conditions holds. In fact, in [S] Bedford and Smillie 
show that if / is hyperbolic, then 51 and f2 are both equal to J union finitely many 
attracting periodic orbits. Thus for hyperbolic polynomial diffcomorphisms of C 2 , 
the basic sets are J and the attracting periodic orbits. 
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3. Interval arithmetic 



In order to genuinely prove dynamical properties, we use in Hypatia a method 
of controlling round-off error in the computations, called interval arithmetic with 
directed rounding (IA). This method was recommended by Warwick Tucker, who 
used it in his recent computer proof that the Lorenz differential equation has the 
conjectured geometry (j^7j). 

In fact we use IA not only to control error, but we take advantage of the structure 
of this method in our algorithms and implementation. We thus give a very brief 
description of IA below, and refer the interested reader to [201 H3 OH ■ 

On a computer, we cannot work with real numbers, rather we work over the 
finite space F of numbers representable by binary floating point numbers no longer 
than a certain length. For example, since the number 0.1 is not a dyadic rational, 
it has an infinite binary expansion. The computer cannot encode this number 
exactly. Instead, the basic objects of arithmetic arc not real numbers, but closed 
intervals, [a, b], with end points in some fixed field K. Denote this space of intervals 
by IK. The operation of addition of two intervals [a, 6], [c,d\ G IK is defined by: 
[a, b]+[c, d] = [a + c, b + d] . Hence if x G [a, b] and y G [c, d], then x+y G [a, b]+[c, d]. 

The other operations are defined analogously, for example: 

[a, b] — [c, d] = [a — d, b — c] , and 

[a, b] x [c, d] — [min(ac, ad, be, bd), max(ac, ad, be, bd)] 

However, an arithmetical operation on two numbers in F may not have a result 
in F. Thus to implement rigorous IA we use the idea of directed rounding to round 
outward the result of any operation. For example, 

[a, b] + [c, d] = [la + c| , T& + dt] , 

where jxj is the largest number in F strictly less than x [i.e., x rounded down), 
and ]x] is the smallest number in F strictly greater than x [i.e., x rounded up). 
For any i£l, let Hull(x) be the smallest interval in IF which contains x. Thus, if 
x G F, then Hull(» = [x, x\. If x G R \ F, then Hull(s) = [[x[ , ]x\}. 

Thus, if the user is interested in a computation involving real numbers, then 
IA with directed rounding performs the computation using intervals in IF which 
contain those real numbers, and gives the answer as an interval in IF which contains 
the real answer. In higher dimensions, IA operations can be carried out component- 
wise, on interval vectors. 

One must think carefully about how to use IA in each arithmetical calculation. 
For example, it can create problems by propagating increasingly large error bounds. 
Iterating a polynomial diffeomorphism like H ac on an interval vector which is not 
very close to an attracting period cycle will give a tremendously large interval vector 
after only a few iterates. That is, suppose B = [a, b] x [c, d] is an interval vector in 
C, and one attempts to compute a box containing P^°(B), for P c (z) — z 2 + c, by: 
for j from 1 to 10 do 
' B = P C (B) 
then the box B will likely grow so large that its defining bounds become machine 
oo, i.e., the largest floating point in F. Similarly, one would also never want to try 
to compute Db„H o • • ■ o Db 1 H o Db H(u), for a vector u G C 2 , since the entries 
would blow up (see Algorithm l5.il StcpQ. 
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Our construction involving boxes in C 2 as the basic numerical objects is designed 
to be efficiently manipulated with IA. For all of our rigorous computations, we use 
IA routines provided by the PROFIL/BIAS package, available at |35j . 

4. Characterizing box hyperbolicity 

In this section we define box hyperbolicity for a box chain model T, in Defini- 
tion l4.il and show that if T is box hyperbolic, then / satisfies the standard definition 
of hyperbolicity, proving Theorem 11.21 Further, we show in Proposition 14. 1 fl that 
box hyperbolicity is equivalent to a simple condition in linear algebra. Throughout 
this section, assume / is a polynomial diffcomorphism of C 2 , with d(f) > 1, and 
let r = (V, £) be a box chain model of / on J. 

Definition 4.1. Suppose for each box Bk in V, we have some nondegenerate in- 
definite Hermitian form, Qk : C 2 — > R. Define Cjt :— {w: Qfe(w) > 0}, as the 
unstable cones, and define their complements as the stable cones: C| := C \ Cj£ = 
{v:Q fc (v)<0}. 

We say that T is box hyperbolic if Df (Df -1 ) preserves and expands the unstable 
(stable) cones, with respect to {Qk}, i.e., for every edge (k,j) £ £(r), and every 
zeB k : 

(1) if w e CI, then D z f(w) E Cf and Q j (D z f(w)) > Q fc (w); 

(2) if v e C], then [D z f]- X (v) e C s k and -Q k {[D z f]- X (v)) > -Q,(v). 

In fact, a given cone Cjjf determines an associated Hermitian form Qk up to 
scaling. Finding an appropriate choice of scaling for each Qk is how we determine a 
metric for which H is hyperbolic. To prove Theorem 1 1.21 we first use a partition of 
unity argument to smooth out the forms {Qk\k=i m ^° a continuous field of forms 
{Qz}ze.j (Definition 14. 8|) . then define a riemannian metric induced by {Q z }zeJ 
(Definition 14.10(1 , and show that H is hyperbolic on J in this new riemannian 
metric. 

4.1. Box hyperbolicity implies hyperbolicity. Here our goal is to prove Theo- 
rcm ll.2l that if T is box hyperbolic, then / is hyperbolic on J, as in Definition ^. II 
Part of the proof is very similar to the one dimensional analog, proved in |19j . in 
that we use a partition of unity to smooth out a discrete norm. But before we 
deal with the norm, we verify that box hyperbolicity implies the existence of a 
continuous splitting preserved by the map. 

Lemma 4.2. If F is box hyperbolic, then there exists a splitting of the tangent 
bundle T Z C 2 = E z © E z , for each z in J, which varies continuously with z in J, 
such that f preserves the splitting, i.e., D z f(E^) — Ej z , and D z f(E") = Ej z . 
Further, for each z E B k , E u z C C™, and E s z C Cf.. 

Proof. Recall that Ncwhousc and Palis show that a diffcomorphism / is hyperbolic 
if there is a field of cones C z (not necessarily continuous) which is preserved and 
expanded by Df, such that the complements are expanded by Df" 1 . In their proof 
f|3U|'). they first show that the existence of a cone field preserved by Df implies 
the existence of a continuous splitting preserved by /, with the unstable (stable) 
directions lying inside the unstable (stable) cones. Box hyperbolicity gives a cone 
field preserved by Df. Thus we have cones C z = C^, if z is in box B k (make 
some consistent choice of box B k containing z, for the benefit of points on the 
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boundaries of the closed boxes). Thus by the proof in j20|, we have the existence 
of the continuous splitting preserved by Df. □ 

It is more natural for computer calculations to use the L°° metric on K 2 " = C", 
rather than euclidcan. Thus, throughout the rest of this paper, ||-|| will denote this 
norm, i.e., if x = (x±, . . . , x n ), then 

||ir|| = max{|Re(x fc )| , |Im(a;fc)| : 1 < k < n}. 

If B is a set of points in C 2 , we denote by N(B, S) the (^-neighborhood of the set B 
in this metric. 

When we say box, we mean a ball about a point in this norm. Thus a box is 
also a vector of intervals, so boxes are neighborhoods which are easily manipulated 
with interval arithmetic. 

We will need to measure the angle between pairs of lines in C 2 , like E% and E z , 
or E z and E%. To do so, we view the set of lines through the origin in C 2 as the 
projective space CP 1 = C. Then the spherical metric on CP 1 induces the following 
metric. 

Definition 4.3. If v = \ Vl ] , w = \ Wl ] are vectors in C 2 , define the distance between 
the directions Cv, Cw to be 

, s . -1 (\V\W2 -V 2 Wi\ 

CT (v, w = sin — — 

V ll v ll ll w ll 

Note that ct(v,w) G [0, tt/2], and for any complex numbers a,/3, 

ct(v, w) = cr(av,/?w). 

Think of this metric as measuring the angle between the complex lines. If v is a 
vector and W is a collection of vectors in C 2 , we define u(v, W) = inf{cr(v, w) : w G 
W}. 

In the next two lemmas, we quantify how, in the metric a, the unstable (stable) 
lines from the splitting of Lemma FOl arc strictly inside the unstable (stable) cones. 

Lemma 4.4. Let T be box hyperbolic. Then there exist d\ > and 5\ > such 
that if ze JC\ N(B k ,5i), then a{E^,C^) > di and a{E s zl C^) > d x . 

Proof. First, note that by compactness of J and the fact that the line fields are 
contained in the interior of the cones, there exists a do > such that 

d < min{cr(E2, C s k ) : z E J n B k } and d Q < min{a(E s z , CI) : z e J n B k }. 

Let d\ = do/2. Next, by compactness of J and continuity of the splitting, there 
exists a Si > such that for any x,z € J with \\x — z\\ < Si, we have a(E^,E^) < di 
and a(E s z ,Ex) < di. 

Now let z e J n N(Bfc) Si). Since z is not necessarily in B k , let m be such that 
z G B mi and a; is a point satisfying x G J n B m n B k and ||x — z\\ < Si. Then 
a(E^,E^) < di and a(E s zl E s x ) < di. Since x G B k we have cj{E l ^C s k ) < d a and 
a{E s x ,Cl) < d . Hence, a(E%, C s k ) > do - di = di, and a(E s zl C^) > d x . D 

Lemma 4.5. Let Y be box hyperbolic. Lf B k ,Bj G V and z G J satisfies z G 
yi(B k ,S x ) and f(z) G ^(B^Si), then a(EJ z ,C^) > d x and a(E s fz ,Cf) > di. 

Proof. This lemma follows directly from Lemma l^H applied to f(z) instead of z. □ 

Before the next step, we need a lemma from |18j . 
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Lemma 4.6 (|18jl. There exists an r\ > so that if Bk,Bj £ V, with z £ INf( J3^ , 77) 
and f{z) £ N(i3j,77), then there is an edge from Bk to Bj in F. 

To prove this lemma, we used the assumption that / was a polynomial mapping 
of degree d > 1, and the fact that by Definition ll.il there is a S > such that there 
is an edge from Bk to Bj if a (5-ncighborhood of f(Bk) intersects Bj. Now, we get: 

Lemma 4.7. Let T be box hyperbolic. Then there is a r > such that for any 
Bk, Bj £ V and any z £ J such that z £ J4(Bk, r) and f(z) £ N(Bj, t), we have 

(1) ifw £ E», then Qj (U,/(w)) > Q fc (w); 

{2) ifv£E s fz ,then~Q k ([D z f]- 1 (v)) > -Qj(v). 

Proof. Among additional requirements given below, let t be less than 77 from 
Lemma 14.61 Then for any z £ J such that z £ 3\f(_Bfc,r) and f{z) £ Ji{Bj,r), 
there is an edge in T from Bk to Bj, i.e., (k,j) £ £. 

Since J is compact, and by Lemmas 14 .41 and 14.51 there exists d 2 > such that: 

d 2 < mm{Qj(D x f(u x )) - Q k {u x ): x £ B k , (k,j) £ £,u x £ E%, \\u x \\ = 1}. 

Let e = ^2/3. By continuity of D x f and the splitting, there is a r < r\ so that 
for any x,z £ J with \\x — z|| < r, z £ N(Bk,r), and f(z) £ "N(Bj,t), we have 

IQfc(uz) -Qk(Ux)\ < e, 
\Q J (D x f(u x ))-Q J (D x f(u z ))\ < c, and 
|Q J (^/(u 2 ))-Q J (C 2 /(u z ))| < e. 

Then Qj(D z f{u z )) - Q fe (u 2 ) > d 2 - 3e > 0. 

Now since Q(w) is a Hermitian form, Q(a~w) = |a| Q(w) for any a £ C. Thus 
by linearity of Df, the above result for u z implies the same result for any w C E™. 
Hence we have Condition (1). 

The proof of (2) is analogous. Let c?3 > satisfy: 

d 3 < min{Qj(s fz ) - Qk(D z r\s fz )): x £ B k , (k,j) € £,s /z £ E} z ,\\s fz \\ = 1}. 

Let e' = ^3/3. Then further restrict r so that for any x, z £ J with \\x — z\\ < r, 
z £ N(B fc ,r), and f(z) £ N(Bj,t), we have 

\Qj(sf z )-Qj(s fx )\ < e', 
\Qk(D x f- 1 ( Sfx ))~Q k (D x f- 1 (s fz ))\ < e', and 
\Q k (D x r 1 (s fz ))-Q k (D z r 1 ( Sfz ))\ < e'. 
Thus (2) follows from Qj{s fz ) - Q k (D z f~ 1 (s fz )) > d 3 - 3e' > 0. D 

Now we use a partition of unity to smooth {Qk}{^ = i on the invariant line fields. 

Definition 4.8. Let T be box hyperbolic. Let r > be as in Lemma l4~7l Define a 
partition of unity on 25 by choosing continuous functions pk '■ C 2 — > [0, 1] for each 
box Bk £ V, such that supp(pfe) C K(Bfc,r) and J2k Pk{ z ) = 1 5 f° r an Y ^ G 23. 
Let zeJ. Then we define Q z : E% U E s z -> R by 

Q,(w)=^p fc (z)Q fc (w). 
fe 

Note that Q 2 (w) is a continuous function of w since Qk is continuous, and 
further a continuous function of z due to the partition of unity. 
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Proposition 4.9. Let T be box hyperbolic. Let {Q z } be given by Definition \4-8\ 
Then for any z G J we have: 

(1) i/w G E%, then Q fz (D z f(w)) > Q z (w); 

(2) ifveEf z: then-Q z ([D z f]- l {v)) > -Q fz (v). 

Proof. Let \i z G E™ be such that \\u z \\ = 1. If wc set 

Qf,z = min{Qj(Dzf(u z )): f(z) G swpp(pj)}, and 

q z = max{Q t (u 2 ):z€ supp(/9 fc )}, 

then by Lemma 14.71 wc know that qj z > q z . Thus wc need only use that the 
partition functions sum to one to get 

Q fz (D z f(u z )) = X)p j (/(«))Q i (I>./(u,)) > X>j(/(*))«3U = ?/,* 

k k 

Hence (1) follows since Df is linear, and for any a G C, Q(aw) = \a\ Q(w). 
Establishing (2) is analogous. Let S/z G P? z be such that ||s/ z || = 1. If we set 

-q s z = mm{-Q k ({D z f}~ 1 (s fz )): z G supp(p fc )}, and 

-Qf,z = max{-Qj(s/ z ): /(z) G supp(ft)}, 

then by Lemma 14.71 we know that — <j| > — g 1 ? „. Thus we need only use that the 
partition functions sum to 1 to get 

-QAlDzf}- 1 (s fz )) = -Y,P^)Qk{[D z fY 1 (s /2 )) > -£>(*)£ = -<£ 



k 



> -iiz = -X>(/(*m, > -E^(/^))^( s /») = -Qf*w- 



Definition 4.10. Let T be box hyperbolic. Let z G J. Wc define the norm 
on T Z C 2 using Q 2 and the spaces E z , E z as a basis, i.e., for w G T Z C 2 , 

||w|| p . 2 =max(|Q 2 (P^(w))| 1/2 ,|Q 2 (P^(w))| 1/2 ), 

where P? z denotes the projection onto E z with P z as its Null space. 



□ 



This is a continuous norm for z G J. Robinson (|36|) notes in his construction of 
an adapted metric for hyperbolic sets that the maximum of two norms on subspaces 
defines a norm, which is very similar to the above. 

Finally, we establish Theorem 11.21 by showing that for a box hyperbolic T, we 
have / is hyperbolic on J with respect to the norm ||-|| z on T Z C 2 , for z G J: 

Proof of Theorem \l.'t& Suppose T is box hyperbolic. We want to show / is hyper- 
bolic over J, as in Definition 12.11 i.e., there is a constant A > 1, and for each z 
in J there is a continuous splitting of the tangent bundle T Z C 2 = E z © E z , and a 
continuous norm ||-|| such that: 

(1) / preserves the splitting, ie., D z f{Et) = E s fz , and D z f{Ef) = E u fz , and 

(2) £)/ (Df^ 1 ) expands on E z (E z ) uniformly, i.e., 

(a) if w G El then \\D z f(w)\\ pJz > A ||w||^, and 
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(b) if w G E% then ||^/- 1 (w)|| p/ _ 1(2) > A ||w|| P)Z . 

We have (1) by Lemma [4.21 Let z G J and ||-|| be given by Definition 14. 1UI 
We show that (2) follows easily from Proposition l4.9l 

First suppose w C E z . Then D z f(w) C E'%. Hence ||w|| = Q z (w) and 
\\ D zf(w)\\lf z = Qfz(D z f(w)). Thus Condition 1. of Proposition Ol implies that 
\\D z fW\\ pJz > \\w\\ PfZ . 

Now consider w C E s z . Then D z f(w) C E s fz . Hence ||w||^ z = -Q z (w) and 
\\D z f(w)f fz = -Q fz (D z f(w)). Then Condition (2) of PropositionOl applied to 
v= (D z f Jl (w)) implies ||D z /- 1 (w)||^ / _ 1(z) > ||w||^. 

Finally, by compactness of J the strict inequalities imply the existence of some 
constant A > 1, proving (2). □ 

4.2. Using linear algebra to characterize box hyperbolicity. First recall 
that a Hermitian form Q : C 2 — ► K is associated to a Hermitian matrix A, so that 
Q(w) = w'iw. Note that if (k,j) G £ is any edge in the graph T, then for any 
z G Bk, (Qj o D z f) is also a Hermitian form, given by 

(3) Qi(D,f(w)) = ^*((D z frA 1 (D z f))w. 

Proposition 4.11. Suppose {Qk} are Hermitian forms with Cj* = {w: Qfc( w ) > 
0} and Cf. = {v: Qfe(v) < 0}, for each box Bk in V. Then T is box hyperbolic 
(using {Qk}) iff for every Bk G V, every z G Bk, and every edge (k,j) G £, we 
have ((Qj o D z f) — Qk) is positive definite. 

Proof. (<=) We begin with the reverse implication. Let z G Bk and Bj be a box 
such that (k,j) G £. Then Qj(D z f)(w) > Q fe (w), for all z G B k and all w G C 2 . 

First consider the unstable cones. Suppose w G C%, so by definition < Qfe(w). 
But then by hypothesis, we get 

0<Qfc(w)<Q,-(I>,/(w)). 

Thus D z f(w) G C", so the unstable cones are preserved by D z f, and we have 
established Condition 1 of box hyperbolicity. 

Next consider the stable cones. First, we show that stable cone preservation 
follows from unstable cone preservation, since they are complementary. Indeed, 
above we showed that Df preserves the unstable cones, i.e., D z f(C^) C C". Hence, 
CI C [D z f]~ l (Cf). But by definition, C s = C 2 \C U . Thus C s k D [D z f\- l {G*) and 
so the stable cones are preserved by Df~ x . 

Now let v G C|, so that 

< -Qj(v) = -(Qj o D z f) ([£>*/]-») . 
Then since we have stable cone preservation under Df -1 , we also know that 

0<-Q fc ([^ z /]- 1 (v)). 

Combining this with the negative of the hypothesis establishes Condition 2 of 
box hyperbolicity, i.e., 

-Qk([D z f]-\v)) > -Qj(v). 
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(=>■) Now we prove the forward implication. Suppose T is box hyperbolic, i.e., 
we have Conditions 1 and 2 of Definition 14.11 Let (k,j) G £, and z G Bk- We 
consider w in each of three regions to show ((Qj o D z f) — Qk) is positive definite. 

Case 1: Suppose w G Cjf. Then by definition < Qfe(w). 

Since box hyperbolicity implies the unstable cones are preserved by Df, we have 
that D z f(w) G Cf, so < Qj(D z f(w)). 

Then Condition 1 of box hyperbolicity gives us 

Qj(D z f(w)) > Q fc (w), 

hence ((Qj o D z f) — Qk) is positive on G u . 

Case 2: Suppose w G [Dgf]' 1 ^), i.e., v = D z f(w) f= Cj. Then by definition 
Qj(D x f(w))<0. 

Now by stable cone preservation, we know w E C|, hence Qfe(w) < 0. 

Condition 2 of box hyperbolicity says that 

-Q,(v) < -Qk([D z f\-Hv)) 

for all vectors u E C s , hence it applies to v = D z f(w). Thus we get 

-Qj(D z f(w)) < -Q fe (w), 

and negating yields 

Q,(D z f(w)) > Q k (w), 

so ((Qj o D z /) - Q k ) is positive on [D z f]- 1 (G S ). 

Case 3: For the remaining w, we have w ^ C^ and w ^ [D z f]~ l (C s A. Then 
Qfe(w) < and Qj(D z f(w)) > 0. Hence, 

Qj(D z f(w)) > > Q fe (w). 

Thus we easily get Q j (D z f(vr)) - Qfc(w) > 0. D 

5. Verifying box hyperbolicity: the Axis Metric Algorithm 

In this section, we explain in detail the Axis Metric Algorithm for testing box 
hyperbolicity of a box chain model T of a Henon mapping, H, by attempting to 
construct a cone field and norm for which the map is hyperbolic. 

But first, before we can test box hyperbolicity, we must start with a T which 
seems to model H reasonably well. Thus we now summarize how we use the box 
chain construction of ^Hl to obtain separate strongly connected graphs modeling 
J and any other invariant sets of recurrent dynamics, for example, sink cycles 
(attracting periodic orbits). Recall from Section l2~3l that if H is indeed hyperbolic, 
then the only recurrent dynamics are J and a finite number of attracting periodic 
orbits. The construction is an iterative process. We begin by defining a large box 
25o = [— R, R] 4 in K 4 = C 2 , such that all the recurrent dynamics of the map is 
contained in So (in |18j we give a simple formula for R in terms of the parameters 
a,c). Then for some n > 1, we place a 2™ x 2" x 2™ x 2™ grid of boxes on "B . The 
construction then builds strongly connected graphs Ti, T 2 , . . . , Tn, each consisting 
of a subcollcction of these grid boxes, and such that J is covered by the boxes of Ti , 
and each sink cycle is covered by the boxes of some Tk- Then T — Ti is a box chain 
model of J, as in Dcfinition ll.il Using smaller boxes in the construction produces 
a more accurate box chain model. 

To prove hyperbolicity, we need each sink cycle in a different model from J, i.e., 
in some 1^ for k ^ 1. If there seems to be a sink cycle together with J in Ti , then we 
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subdivide and repeat the above process. That is, place a grid of boxes inside of each 
box of Ti, and use these smaller boxes to obtain a refinement, I\i, I\2, ■ ■ ■ , Fi,M, 
such that r = I\i contains J. If in this refinement, the sink cycle is in some i\i 
for k > 1, then we can stop and study I\ Otherwise, repeat the subdivision process, 
until computational resources are exhausted, or a T containing only J is produced. 
We can check our accuracy at each level in the iterative process by producing 
pictures of the current T's boxes intersected with an unstable manifold of a saddle 
periodic point. As discussed in Section r2.2l we can parametrize an unstable manifold 
by a plane, then to determine the coloration of a pixel, we check whether the pixel 
intersects some boxes of T. Since the picture is a parametrization of a manifold 
which does not line up with the axes in C 2 , a pixel may hit more than one box, 
and in more than one strongly connected component Tfc. The user may also decide 
to lighten the pixels which are heuristically found to be in K + , to check visually 
how close the model is to J. For example, for the Hcnon mapping H ac , with 
(c,a) = {— .3, .1), Figure [3] shows a parameterized unstable manifold intersected 
with the boxes in models ri,r2,r3, with Ti containing J, T2 containing the fixed 
sink, and T3 containing pseudo-recurrence but no true recurrence (thus T3 would 
be eliminated for smaller box size). In this figure each Tfc is shaded differently, and 
pixels heuristically found to be in K + are lightened. 

Now assume we have obtained a model T which appears to contain J but no other 
recurrent dynamics. We outline below the Axis Metric Algorithm for verifying box 
hyperbolicity on T, then describe each step in detail. 

Algorithm 5.1 (Axis Metric Algorithm). (1) Define an approximately invari- 

ant splitting; specifically, a pair of "unstable" and "stable" vectors, Ufc,Sfc, 
for each box, B k e V(T). 

(2) Build an "unstable" metric which is (approximately) uniformly expanded 
by DH on the set of unstable directions, {Cu^}, and a "stable" metric 
which is (approximately) uniformly expanded by DH -1 on the set of stable 
directions, {Cs^}. 

(3) Use the directions and their metrics to define the cones as Hermitian forms. 

(4) Finally, check whether H preserves the cone field, and whether with respect 
to the Hermitian forms, DH (DH~ l ) is expanding on the unstable (stable) 
cones. 

In this algorithm, we will need to move around on the graph T using the following: 

Definition 5.2. A directed graph Tq is an arborescence if there is a root vertex vq 
so that for any other vertex v, there is a unique simple path from vq to v. Such a 
graph is a tree, and must have exactly one incoming edge for each vertex v ^ wo- 
lf r is strongly connected, then for each vertex vq in T, there is a minimum span- 
ning tree Tq with root vertex vq which is an arborescence (simply perform a depth 
first or breadth first search from vq). We call such a To a spanning arborescence. 
(See 2H| for a discussion of implementation of basic graph theory). 

Step 1 (Setting stable and unstable directions). Recall from Section l2~2l that when 
\a\ y^ 1, except on the curve of equation 4c = (1 + a) 2 , the map H a ^ c has at least 
one saddle fixed point, p. If H has no saddle fixed point, it should be possible to 
instead use a saddle periodic point of period greater than one, but this was not 
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Figure 3. Shown above is the unstable manifold parameterization 
for a saddle fixed point of -ff a ,c, with c = — .3, a = .1, intersected 
with three box models ri,r2,r3, containing pseudo-recurrent dy- 
namics of H. The two islands in the large interior on the left are 
in T2, the model of the fixed sink. The right hand band is Ti, 
containing J. The lighter band on Ti is approximately contained 
in K + . The darkest spots overlapping Ti show T3, which con- 
tains pseudo-recurrence but no true recurrence. Thus T^ would 
be eliminated for smaller box size. Here boxes are of side length 
(2 x 1.43)/2 7 = 0.022, and Ti contains 32,000 boxes. This T x is 
box hyperbolic. 
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necessary for any of the maps we were interested in studying. The Henon mapping 
(x, y) — ► (x 2 + c — ay,x) has two fixed points. Note first that fixed points of H 
must be on the diagonal, i.e., x = y. Then the two fixed points are: 



l -([l + a)±^{l + a) 2 +Ac) 



x+ 
± 2 

By substituting the IA Hull of a and c into the above formula, and performing 
operations in IA, we compute x± as complex intervals (interval vectors in R 2 ) 
containing the actual fixed points. Next, the eigenvalues of D X H for the fixed point 
(x, x) are: 

A, [i — x ± v x 2 — a. 

Thus we test to see whether each interval vector x± is a saddle, by computing two 
real intervals containing the moduli of the eigenvalues, and then testing whether 
one interval is entirely greater than one, while the other is entirely less than one. 
So now assume H has a saddle fixed point p = (x,x). Bedford and Smillic have 
shown J contains all saddle periodic points, hence p £ J, thus there is a box, Bq, 
in r, containing p. Let Uo (so) be the eigenvector of unit length corresponding to 
the unstable (stable) eigenvalue. These are natural unstable and stable directions 
in the box Bq. Since we will just use the Uk's and Sk's as axes for cones, we need 
only know them approximately, so here interval arithmetic is not needed. 

Now let Tq be a spanning arborescence of T, with root vertex Bq. To define 
unstable directions in each box B k , sucessively push Uo across the edges of Tq by 
DH . To be precise, fix a point Zk in B k , say the center point of the box, and for 
each edge (k,j) £ £(ro), starting with fc = 0, define 

D Zk H(u k ) 



\D Zk H(u k ) 



As noted above, we need only an approximation to the u^'s, so interval arith- 
metic is not used here. In fact, were we to attempt to use IA, we may encounter a 
computational difficulty. This is due to the recursive definition of the u^'s. Suppose 
we desired to compute "more" of the potential unstable directions in each box, by 
starting with an interval vector Uo in C 2 guaranteed to contain the unstable eigen- 
vector of D p H , and then pushing this interval vector across Tq to define interval 
vectors Uk in each box B k . More precisely, a box B is an interval vector in C 2 , so 
B = (X, Y), with X, Y complex intervals (interval vectors in M 2 = C). Then 



D p H = 



2p —a 
1 



and DuH = 



2X -Hull(o) 

Hull(l) Hull(O) 



so DbH has entries complex intervals. Then DB k H(Uk) is an interval vector. 
Suppose (1,2,..., n) is a path in Tq (this should be (fci, . . . , k n ), but we wish to 
avoid too many nested subscripts). Then U n = Ds n H o ■ ■ ■ o Db 1 H(XJ 1 ). However, 
since the B k s arc near J, due to the dynamical properties of H on J, this type 
of iteration will (in experiment, very quickly) lead to intervals U n to huge to be 
useful. 

Define the stable directions similarly, keeping in mind that stable cones should 
be expanded and preserved by DH^ 1 . The transpose of a graph T, T T , is the graph 
formed by reversing the edge directions of I\ Thus to define Sfc, we use a spanning 
tree Tq of T T with Bq as root vertex (note that Tq ^ Tq), and push Sq across 
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successive edges of T by DH~ l . Specifically, if (j, k) e £(r ), then 

D Zk H-\s 3 ) 
k \\D Zk H-i( Sj W 
where z& € B^ as in the definition of Uj. . 

Proposition 5.3. Let {u^}, {s^} be unstable and stable directions for each box Bk 
ofT, defined as above. 

(1) Let C% be any box constant cone field preserved by D Z H, for each z £ B^. 
Then for each Bk, we must have \ik C C%. 

(2) Let Cf. be any box constant cone field preserved by D z H~ l , for each z £ -Bfc. 
Then for each Bk, we must have S& C C|. 

Note that for edges (k,m) in T but not in the spanning arborescence Tq, DH 
does not map U& into Cu m . It is helpful in establishing invariancc of the cone field 
if DH(\ik) is close to Cu m . and in addition, if D Z H does not vary greatly as z 
varies within one box. 

Thus before performing the next step, we take some measurements on the sepa- 
ration of the stable and unstable directions in each box, to get an idea of whether 
it might be possible to prove box hyperbolicity using these directions, and store 
(for later scrutiny) which boxes might be obstructions. In order to measure the 
difference between directions, we view a direction in C 2 as a complex line in C 2 , 
and thus use the spherical metric, a. Then, for each box Bk in V(T), let 

Udiam[fc] = dianv{-D^-ff(uj): (j,k) 6 £(T)}, and 

Sdiam[fc] = diam a {[D Zk H]- 1 (sj) : (k,j)e£{T)}. 

We do not measure the variation within one box, i.e., between D Z H and D Zk H 
for different z in box Bk, since it seems that would be much smaller than among 
images from different boxes. 

Proposition !.^. 31 suggests that a clear separation between Udiam[fc] and Sdiam[fc] 
is needed in order for a computer program to verify cone preservation under DH, 
thus we will not confidently progress to the next step unless we have 

c-(ufc,Sfc) — (Udiam[fc] + Sdiam[fc]) > 

in each box Bk- Finally, wc record in which boxes the either Udiam or Sdiam is 
large, or a — (Udiam + Sdiam) is negative. 

Step 2 (Building a metric on the directions). Consider the unstable directions 
{C\ik}k =1 . As discussed above, DH does not quite preserve these directions, so 
first we take that into account. Let P s u be the projection onto Cu with Cs as its 
Null space. Given the vectors u = (ui,U2) and s = (si,S2) in C , this projection 

is: 



P. 1 



ViS 2 -v 2 si 



UlS 2 - U 2 Sl 

Then for each edge (fc, j) <E £(r), we have P s Uj o D Zk H maps Cu,t to Cuj. If the 
unstable directions are a good approximation to an invariant unstable line field, 
then P s u o DH is close to DH on these unstable directions. 

In |19| , we describe a method for proving hyperbolicity of polynomial maps of 
one complex variable, by building a metric for which the map is expanding by some 
L > 1 on a neighborhood of the Julia set. The neighborhood of J is a collection 
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of boxes in C, and the metric in each box is defined by a constant times the L°° 
metric in R 2 = C. Say tfk 1 1 - 1 1 is the metric on box Bk, then the constants <pk are 
called metric handicaps. 

In two variables, we use this same algorithm twice: once for the unstable direc- 
tions and once for the stable directions. For example, for the unstable direction 
field, we will attempt to build a metric for which P s u o DH is box expansive by 
some amount L > 1 . This metric will be defined in each unstable direction Ciik by 
a constant times our L°° metric in R 4 = C 2 , say tp% ||-||. Following the algorithm in 
|19j . we call the tp% unstable metric handicaps, and define ipfi = 1 in box B$, then 
want to build handicaps satisfying 

^\\P£oD Zh H{v)\\>Lipt\\v\\, 

for each edge (k,j) G £(r), and each v G Ciik. But then v = auk for some a G C, 
so since P s u o DH is linear and ||ufc|| = 1, the above equation is equivalent to 

( A \ ,u ^ L tp% 

4 Vj > jr^Tj „ v , TTT- 

If we set £k,j = £/||-Ps" : ° Dz k H(uk)\\, then we can use Algorithm 4.4 of ^H] to 
attempt to find metric handicaps satisfying ip™ > Cfe.iVfe) hence Equation 21 on 
each edge (k, j) G I\ This is of course not always possible, but the intuition is that 
it should be possible if the box model is sufficiently small in the right places. If it 
is not possible, then the user can try a smaller L, or start over with smaller boxes. 

We use interval arithmetic in Equation 0] as in ^Hj • The <p% 's will be chosen in 
F (machine knowable numbers, note tp^ = 1 G F). Then replace each number in 
the left hand side with it's Hull (the smallest interval in IF containing it). Then 
perform the operations in IA, and test whether the upper endpoint of the resulting 
interval is less than tp™. 

To attempt to define stable metric handicaps, use the method analogous to 
that for the unstable metric handicaps. That is, try to find an M > 1 and build 
handicaps {<p s k } so that 



(5) V* < 



\P&oD Zh H(a k ) 



for each edge (k,j) G £(r). Then the stable directions are box-contracted by 
P„ o DH. Again, if this step fails to produce a contracted metric, then the user can 
try a larger M, or start over with smaller boxes. 

As we will see below, the ratio of stable to unstable metric handicaps in each 
box, (pf. to iffr, determines the width of the cones. Hence it is necessary to find 
values for L and M which yield comparable metrics. Lyapunov exponents give us 
some intuition. For a polynomial automorphism of C 2 , there are two Lyapunov 
exponents, A , which measure expansion and contraction of tangent vectors. Ac- 
cording to Q, if / is a polynomial diffeomorphism of C 2 with d = d(f) > 1, then 
A + > logd, A~ < — logd, and 

(6) A + + A _ =log(det£>/). 

Note that for Hcnon mappings, det(£>/) = a. Thus A + + A~ = log(a). Since 
A+ > logd, we have the inequality: A - < log(a) — logd. In the case \a\ < 1, we 
have log(a) < 0, hence the inequality for A~ is stronger than the inequality for A + . 
Thus in general we expect stronger contraction than expansion of tangent vectors 
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under Henon mappings. Equation El implies that a good rule of thumb for choosing 
L and M is LM = \a\. In practice, for any T we tested, the algorithm for setting 
the stable metric handicaps given any M always completed in much less time than 
the algorithm for setting unstable metric handicaps given an L (perhaps due to 
the strong contraction). Thus after experimenting with various options, we have 
adopted the following strategy. First find the smallest working M using simple 
bisection (0 < M < 1), then test expansion on the unstable directions using a value 
of L near \a\ jM. 

In experiment we have observed that finding a working L > 1 and M < 1 is 
almost always possible when a box model T has been found for J which does not 
contain any sinks of H. Rather, the difficult step is the next one: checking whether 
the cone field defined by these metrics is preserved and expanded by DH. 

Step 3 (Defining a cone field). If Step successfully constructed expanded and 
contracted metrics on the unstable and stable directions, respectively, then the 
metrics and directions always define cones in each box, as follows. 

In each box, Bk, define the unstable cone, Cg, so that a vector w is in the 
unstable cone if it is closer to u^ than Sfe, relative to the unstable and stable 
metrics. That is, w £ C% iff 



^Ik*(w)|| >v*ra(w)i 



Then the stable cones are just the complements, Cf. := C 2 \ C£ 



We define the Hermitian form Qk : C — > R, by 

Q k (w) = (^||P s ^(w)||) 2 - (^||P^(w)||) 2 . 

Thus the unstable cone, C£ , is simply the set of vectors for which Qk is nonnegative, 
and the stable cone, Cf., is the set of vectors for which the form is negative. 

Wc can construct a Hermitian matrix, A k , which encodes the information of Qk, 
following standard linear algebra as in |16j . A Hermitian form Q defines a sesquilin- 
ear form j:C 2 xC 2 ^ M, such that <?(w, w) = <2(w), where we can recover g using: 

1 1 i i 

g(v, w) = -Q(v + w) - -Q(v - w) + -Q(v + iw) - -Q(v - iw). 

A sesquilinear form g can be represented by a matrix A so that g(y, w) = w* Av, 
with a m . n = g(e n ,e m ) for an ordered basis {ei,e 2 }, like {(1,0), (0, 1)}. Now g 
Hermitian implies that A is also Hermitian, and the range of Q is R. Thus, <5(w) = 
w'Aw, where a m . n = |Q(e n + e m ) - \Q{e n - e m ) + \Q(e n + ie m ) - \Q(e n - ie m ). 

We easily calculate that for u = {u\, U2), s = (si, S2), if wc set 

hi = (¥>" \s 2 \ \\u\\) 2 - (^ \u 2 \ || S ||) 2 , b 22 = {<? \ Sl \ \\u\\f - (<p s |«i| ||.s||) 2 , 

hi =4 W 1 ||w||) (|s2 — s\\ -|«2 + si| +i|is 2 -si| -i\is 2 + si\) 

- (tp s \\s\\) 2 (\u 2 - ui\ -\u2 + ui\ +i\iu2—ui\ —i\iu2+ui\) 

and 621 = h 2 , then a m , n = b m . n /\u 1 s 2 - u 2 Si| ■ 

In implementation, we use the above formulas to calculate, for each box B^ in 
r, an interval valued matrix A^, representing the form Qk defining the cone C^. 
Note that since A/, is Hermitian, the main diagonal entries are real intervals, and 
the other entries arc (complex conjugate) complex intervals. As noted above, the 
vectors Uk and Sk and the handicaps <p)f , ipf. are chosen to be machine knowable 
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numbers. Thus before the formulas defining Ak are computed, these terms are 
converted to their interval Hulls, then the arithmetic operations are performed in 
IA to obtain Ak with interval entries (of length greater than zero) . 

Remark 5.4. Note that the ratio of the metric handicaps determines the angle 
width of the cones. Thus if ip^ and (p k are several orders of magnitude different 
then the cones will be very thin, even if the unstable and stable directions are far 
apart, thus the cones will be difficult for the computer to work with. This is why 
it is necessary in Step to find values of L and M which yield a comparable pair 
of metric handicaps in each box. 

Step 4 (Checking whether cones are preserved and expanded). For the last step 
of testing box hyperbolicity, we need to test whether DH (DH -1 ) expands the 
unstable (stable) cones, with respect to {Qk}- For this step we simply use Propo- 
sition 14.111 in which we showed that in order to get preservation and expansion 
of the unstable cones, and contraction of the stable cones, we need precisely that 
((QjoD z H) — Qk) is positive definite for every edge (k, j) S £(r), and every z 6 Bk- 
Thus in this step, we simply compute this form defined for each edge in the graph, 
and test whether it is positive definite. 

In Step 01 for each box Bk in T, we computed an interval matrix Ak representing 
Qk, in that Qfc(w) = w*^lfcw. Thus by Equation|3 the interval matrix representing 
(Qj o D Z H) - Q k is 

T kj = ((D Bk Hy A, D Bk H) - Ak. 

Using the formula for the interval matrix Ds k H from Step ^ it is straightforward 
to compute the interval matrix Tk.j. Now we need to check whether this matrix is 
positive definite. But since T is Hermitian, the trace and determinant are real, and 
T is positive definite iff its trace and determinant are positive (see QB]). The trace 
and determinant of T are real intervals, so we simply compute them with IA and 
check whether their lower endpoints are positive. 

If the above test succeeds (positive definite for each edge), then the model T is 
box hyperbolic, hence H is hyperbolic. If not, then we may record boxes which are 
obstructions, that is, boxes Bk, Bj for which Tkj fails to be positive definite. 

This is the stopping point of the Axis Metric Algorithm for testing box hyper- 
bolicity for a given T. If box hyperbolicity fails, the user may refine T by choosing 
to subdivide either all the boxes, or some subset of the boxes which seem to be 
obstructing the hyperbolicity testing (for example, boxes marked in Steps ^ or 0] 
above), then test the new T with the Axis Metric Algorithm. Two of the Hcnon 
mappings of Theorem 11.31 H ac with (c,a) — (0, — .22) and (1.5, .5), were proven 
hyperbolic for a model constructed by straight subdivision to a certain box size, 
then by subdividing twice only boxes marked in Step 01 (see next section) . 

6. Results of running Hypatia on Henon mappings 

Running Hypatia for a Hcnon mapping is not quite as simple as inputting the 
parameters a, c and awaiting the results, since the user must make decisions as to 
how to build the best T for testing with the Axis Metric Algorithm. In this section, 
we describe the specific process we followed and results obtained for the mappings 
of Theorem IP 
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FIGURE 4. A FractalAsm picture of W™ (~l K + for H a<c , with c 
—3, a = —.25. This set appears to be a real Cantor set. 



Theorem 11.31 follows from Theorem 11.21 that box hyperbolicity of some T, a 
box chain model of J, implies hyperbolicity of H on J, and from the fact that for 
each mapping mentioned in the theorem, using our computer program Hypatia, we 
constructed a box chain model of J and verified box hyperbolicity with the Axis 
Metric Algorithm. Below, we discuss the process for each of these mappings in 
increasing order of the computational difficulty of proving hyperbolicity. 

The quickest map to be proven hyperbolic was H a>c , with c = —.3, a = .1. We 
simply used a box chain model of J with boxes selected from an evenly subdivided 
(2 7 ) 4 grid on T> = [-1.43, 1.43] 4 . Figure shows the box hyperbolic I\ This is a 
map seemingly in the main cardioid, with recurrent dynamics J and a fixed sink. 

For Ha tC , with c = —3, a = —.25, we also proved hyperbolic relatively quickly. 
The box chain model r of J from an evenly subdivided (2 7 ) 4 grid on Bq = 
[—2.57, 2.57] 4 is box hyperbolic. This mapping appears to be a real horseshoe 
(i.e., a horseshoe contained in R 2 ). Figure 0]is a FractalAsm picture of the Julia 
set. This kind of picture is really the most useful for a Cantor set. 

We proved the map H ac , with c = 0, a = —.22, is hyperbolic by starting with a 
model of J from an evenly subdivided (2 7 ) 4 grid on 23 = [-R, -R] 4 , for R = 1.32, but 
then additionally performing three hyperbolicity tests, and each time subdividing 
only boxes in which the cone check of Step 0] of the Axis Metric Algorithm failed 
(to end with boxes of size ranging from 2R/2 7 to 2R/2 10 ). This map also seems to 
be in the main cardiod. The picture of the Julia set is similar to Figure 03 

Using nearly the same method as the previous mapping, we proved H aet with 
c = 1.5, a = .5, is hyperbolic. Here, we started with the even (2 8 ) 4 grid on 
25o = [-R, -R] 4 , for R = 2.286, then twice subdivided only boxes in which the 
cone check (Step 0] of the Axis Metric Algorithm) failed (yielding boxes of size 
2R/2 8 to 2R/2 la ). The resulting box chain model is box hyperbolic. FractalAsm 
pictures (see Figure EJ suggest this map is a complex horseshoe, with Julia set not 
contained in K 2 . 

Table ^ contains more data for all of the mappings discussed in this section. In 
the table, T denotes the box chain model of J, for the map H a>c (x,y) = (x 2 + 
c — ay,x). R is the bound such that the boxes are contained in CE>o = [ — R, R] 4 C 
R 4 = C 2 . The box grid depth for a box is the number n such that the box is of size 
2R/2 n . If a box chain model contains boxes of multiple sizes, then multiple box 
grid depths are listed. 
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FIGURE 5. A FractalAsm picture of W™ (~l K + for H a<c , with c = 
1.5, a = .5. This set appears to be a complex Cantor set. The 
origin of the axes is the saddle fixed point p. 
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